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ABSTRACT 

Aims. In a previous paper we have investigated the molecular environment towards the eastern border of the SNR G18.8-h0.3. 
Continuing with the study of the surroundings of this SNR, in this work we focus on its southern border, which in the radio continuum 
emission shows a very peculiar morphology with a corrugated corner and a very flattened southern flank. 

Methods. We observed two regions towards the south of SNR G18.8-h0.3 using the Atacama Submillimeter Telescope Experiment 
(ASTE) in the ^^CO J=3-2. One of these regions was also surveyed in ^^CO and C^^O J=3-2. The angular and spectral resolution of 
these observations were 22", and 0.11 km s"^ We compared the CO emission to 20 cm radio continuum maps obtain as part of the 
Multi-Array Galactic Plane Imaging Survey (MAGPIS) and 870 pm dust emission extracted from the APEX Telescope Large Area 
Survey of the Galaxy. 

Results. We discovered a molecular feature with a good morphological correspondence with the SNR’s southernmost corner. In 
particular, there are indentations in the radio continuum map that are complemented by protrusions in the molecular CO image, 
strongly suggesting that the SNR shock is interacting with a molecular cloud. Towards this region we found that the ^^CO peak is not 
correlated with the observed ^^CO peaks, which are likely related to a nearby Hii region. Regarding the most flattened border of SNR 
G18.8-h0.3, where an interaction of the SNR with dense material was previously suggested, our ^^CO J=3-2 map show no obvious 
indication that this is occurring. 
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1. Introduction 


The investigation of the surroundings of supernova remnants (SNRs) can provide useful information for a wide range of research 
fields, as the physical and chemical changes induced by the passage of a shock wave, dust production and destruction, the possibility 
of triggering stellar formation, production of high-energy gamma rays, etc. Multiwavelength studies are required to explore these 
effects. 

The SNR G18.8+0.3, with a peculiar shape wi th flattened borders along the east and sout h , has been suggeste d to be interacting 
with dense molecular gas towards such directions (iDubner et al.lll999ll2004tlTian et al.ll2007h . iParon et ^ (l2012h . hereafter Paper I, 
investigated the molecular gas towards the eastern border of the SNR G18.8+0.3. They found a dense molecular clump, close to 
the shock front of the SNR but not in contact with it, containing a complex of Hii regions and massive young stellar objects (YSOs) 
deeply embedded. This region is indic ated in Eig.fjas ‘Eas tern Clump’. Eigure[T] shows the SNR G18.8+0.3 in radio continuum at 
20 cm as extracted fro m the MAGPIS (iHelfand ^ aDl200^ (in blue) and its surroundings as seen at the 870 pm emission extracted 
from the ATLASGAL (iBeuther et al.ll2012l) fin green). Southwards of the SNR, where the molecular cloud appears to be in contact 
with SNR shock front (see Eigure 1 in Paper I), there are a few Hii regions and submillimeter dust compact sources, making this 
region an interesting target to study the interplay between the SNR shock front, surrounding molecular gas, and Hii regions. 

Using molecular lines we investigate two regions towards the south of SNR G18.8+0.3; Reg.l and the subregion Reg.lb, were 
selected to investigate the possible caus es for the peculiar mo rphology of the SNR southernmost ‘corner’, and the molecular ga s 
towards the Hii region GOl8.630+0.309 (lAnderson et al.ll201ll) and the ATLASGAL source 018.626+00.297 (lUrquhart et al.ll2014l) . 
while Reg.2 was chosen to investigate the molecu Ur environment tow ards the most flattened border of the SNR where originally it 
was suggested an interaction with a dense cloud (iDubner et al.lll999[). As in Paper I, her e we assume a distance of 14 + 1 kpc for 
SNR G18.8+0.3. Concerning to the Hii region GOl8.630+0.309. lAnderson et al.l (l2012l) could not resolve the distance ambiguity 
between 1.5 and 14.7 kpc. 
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Fig. 1. Two-color image of the SNR G18.8-I-0.3 and its surroundings. The radio continuum emission at 20 cm (from MAGPIS) is displayed in 
blue (right color bar), and the ATLASGAL emission at 870 fim in green (left color bar). Both color ranges are in Jy beam“^ The surveyed regions 
in the molecular lines using ASTE are indicated. 

Table 1. Surveyed regions using ASTE. 


Region 

Center (R.A., dec.) 

Size 

Reg.l 

18:23:48.0, 

-12:32:14.5 

3'x4' 

Reg. lb 

18:23:47.1, 

-12:33:12.1 

2'x2' 

Reg.2 

18:23:38.0, 

-12:30:04.6 

3'x3' 


2, Observations 

The observations of the molec ular lines were carr ied out between September 7 and 13, 2013 with the 10 m Atacama Submillimeter 
Telescope Experiment (ASTE: lEzawa et all2004l) . We used the CATS345 receiver, which is a two-single band SIS receiver remotely 
tunable in the LO frequency range of 324-372 GHz. The surveyed regions are shown in Eig.[T]and summarized in Table[T] The 
observations were performed in position switching mode and the mapping grid spacing was 20" in all cases. The J=3-2 

line (345.796 GHz) was observed towards Regs. 1,2. In Reg.lb we also observed J=3-2 (329.330 GHz), and J=3-2 

(330.587 GHz). The integration time was 20 sec for the and 50 sec for the other molecular lines in each pointing. 

We used the XE digital spectrometer with bandwidth and spectral resolution set to 128 MHz and 125 kHz, respectively. The 
velocity resolution was 0.11 km s“^ the half-power beam width (HPBW) 22", and the main beam efficiency was T/mb ~ 0.65. The 
spectra were Hanning smoothed to improve the signal-to-noise ratio and only linear or some second order polynomia were used for 
baseline fitting. The data were reduced with NEWSTArQ and the spectra processed using the XSpec software packageEI. 

3. Results and discussion 

3.1. Regions 1 and 1b 

In Reg.l, J=3-2 emission is detected at velocities between 5 and 30 km s“^, very much like what was observed in the eastern 
region (Paper I). The emission integrated over the mentioned velocity range is displayed as contours over a two-color image showing 
the 20 cm and 870 fim continuum images in Eig.O Importantly, the boundaries of the radio continuum emission from the SNR 
and the emission from the molecular cloud show a number of common features, strongly suggesting that the shock front is 
interacting with the cloud. In particular, the most noticeable indentations in the border of the SNR are complemented by protrusions 
in the molecular cloud. In addition, towards the south the contour at 56 K km s“^ of the emission has a peculiar curvature 
apparently in coincidence with the Hii region G018.630-f0.309, suggesting that this object is also associated with the molecular 
structure. 

EigureO displays the J=3-2 emission integrated between 15 and 26 km s“^ towards Reg. lb in white contours. Eor context, 
the 56, and 68 K km s“^ contours of the emission shown in Eig.[2]are also shown in Eig. [3] (yellow dashed contours). The 
emission reveals a molecular feature, composed by two clumps, that partially surrounds the Hit region G018.630-f0.309. One 
of these clumps, the southernmost, exactly overlaps the source ATLASGAL 018.626-f00.297. Curiously the most intense emission 
of the is not correlated with the peaks. Concerning the J=3-2 fine, we only detect emission (a single spectrum) 
over the maximum of the ATLASGAL source, implying that this is the region with highest density. 

In order to discuss the molecular gas in Reg.l and Reg.lb we denote the structure delimited by the 56 K km s“^ contour (the 
dashed yellow contour in Eig.O as the “^^CO-clump”, the feature delimited by a circle of 25" in radius centered at the northern 
maximum in the emission is designated as the “^^CO-clumpN”, and the southern clump delimited by the 16.5 K km 
s“^ contour as the “^^CO-clumpS” (see EigO. We estimated the mass of these features applying a variety of methods, namely: 
by assuming local thermodynamic equilibrium (LTE), from the CO luminosity (in the case of the ^^CO-clump), from the thermal 
dust emission at 870 jim (in the case of the ^^CO-clumpS), and based on the virial theorem. To obtain the LTE mass, we follow 

^ Reduction software based on AIPS developed at NRAO, extended to treat single dish data with a graphical user interface (GUI). 

^ XSpec is a spectral line reduction package for astronomy which has been developed by Per Bergman at Onsala Space Observatory. 
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Fig. 2. Two-color image towards Reg.l with the radio continuum emission at 20 cm (from MAGPIS) displayed in blue (right color bar), and the 
ATLASGAL emission at 870 jim in green (left color bar). Both color ranges are in Jy beam"^ The J=3-2 emission integrated between 5 
and 30 km s~^ is presented in contours, with levels of 32, 42, 56, and 68 K km s“^ The beam of the molecular observations is shown at the top 
right corner. The positions of the Hii region GOl 8.630+0.309 and the compact dust source ATLASGAL 018.626+00.297 are indicated. 



Fig. 3. Two-color image towards Reg. lb with the radio continuum emission at 20 cm (from MAGPIS) displayed in blue, and the ATLASGAL 
emission at 870 fim in green. The color ranges are the same as in Fig.[2| For context, the box of Reg.l and the 56, and 68 K km s~^ contours (in 
dashed yellow line) of the emission presented in Fig.[2|are also shown. Contours of the J=3-2 emission integrated between 15 and 26 
km s~^ are shown in white, with levels of 12.0, 15.0, 16.5, and 20.0 K km s“^ The beam of the molecular observations is shown at the top right 
corner. 


the procedures described in Paper I to estimate the column densities of N(^^CO) and N(^^CO). Tabled presents the derived and 
used parameters to perform that. The column density values represent the total column density as obtained from summation over 
all b eam positions belon ging to each molecular structure described above. To obtain the mass from the CO luminosity we use 
(e.g. iBertsch et al.ll 19931) : N(H 2 ) = 2x 10^^ W(^^CO) to estimate N(H 2 ), where W(^^CO) is the velocity integrated emission 
along the whole ^^CO-clump structure. Finally, once obtained the total N(H 2 ), either from the LTE approximation or from the CO 
luminosity, the mass is derived using M = ju mu ^ N(H 2 ), where Q is the solid angle subtended by the beam size, mu is the 
hydrogen mass, yu, the mean molecular weight, assumed to be 2.8 by taking into account a relative helium abundance of 25 %, and 
d is the distance. 

Assuming that the ^^CO-clumpS is related to the ATLASG AL source 018 . 626+0 0.297, we estimate its mass from the thermal 
dust emission at 870 jum following the procedure presented in ICsengeri et al.l (l2014l) . We use an integrated flux of S^t = 3.31 Jy 
(lUrquhart et al .1120ill) , and assume a typical dust temperature of Td = 15 K. 

Finally, the virial mass of the three molecular structures is derived from the relation: 


Mv 


5 Rclump(Av)2 

8 Gln(2) 


( 1 ) 


where G, Rdump and Av are the gravitational constant, the clump radius (assuming a spherical geometry), and the line width, respec¬ 
tively. Following ISaito et aH (l20Q7l) we define the clump radius as the deconvolved radius calculated from Rdump = 
where S is the area inside the clump. The velocity width Av was derived through a Gaussian fitting to the averaged spectrum 
obtained for the molecular features. The deconvolved radius (Rdump) are 2.95, 1.52, and 1.05 pc, and the Av = 7.05, 3.80, and 5.10 
km s“^ for the ^^CO-clump, ^^CO-clumpN, and ^^CO-clumpS, respectively. All mass results are presented in Tabled 
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Table 2. Used and derived parameters for the LTE mass. 


Parameter 

‘^CO-clump 

‘^CO-clumpN 

^"^CO-clumpS 

Tex (K) 

16.6 

12.5 

11.0 

t‘2 

15 

- 

- 

t 13 

- 

1.3 

3.5 

N('2C0) (cm-2)* 

7.1 xlO^^ 

- 

- 

NC'^CO) (cm-2)t 

- 

1.5 xlO^^ 

3.5 xlO^^ 

N(H2) (cm-2) 

7.1 xl022 

7.5 xl022 

1.7x1023 


Notes. To obtain the N(H 2 ) from the N(^^C O) we assume the c anonical abundance ratio [^^CO/H 2 ] = 10 To obtain the N(H 2 ) from the 
N(^^CO) we assume [^^CO/H 2 ] = 2 x 10"^ (e.g. ISimon et aLl[200lh . 

Table 3. Mass values. 


Clump 

Mlte 

(Mo) 

Mco lum. 

(Mo) 

Mdust 

(Mo) 

Mvir 

(Mo) 

^^CO-clump 

2.7 X 10 ^ 

6.0 X 10 ^ 

- 

3.1 X 10 ^ 

^^CO-clumpN 

3.0 X 10 ^ 

- 

- 

4.5 X 10 ^ 

^^CO-clumpS 

7.0 X 10 ^ 

- 

5.5 x 10^ 

5.7 X 10 ^ 


The only case in which the virial mass differs considerably from the mass obtained by other methods is for the ^^CO-clump. For 
this feature, the virial mass is almost an order of magnitude larger than the LTE and CO luminosity based-estimates of the mass. 
This suggests that this feature is not gravitationally bound. 

We assume as a working hypothesis that the ^^CO-clump is a small molecular feature supported by external pressure. The 
external pressure required for theJ^^CO -clump to remain bound, can be roughly evaluated from the following expression of 
the virial theorem (e.g. iKawamura et al.|[T9^ : 

lU + Q.- ^nRaampPext = 0 (2) 


where 


U = 3/2M 


Av^ 

81n(2)’ 


and 


3GM2 

^I^clump 


(3) 


Using a mean value between Mlte and Mco lum. for the mass, we obtain that it is required a Pextl^B ~ 1.4 x 10^ cm“^ K, with ks the 
Boltzmann c onstant, to keep the ^^CO -clump bound. This value is at least one order of magnitude larger than the external pressures 
calculated in IKawamura et al.l (1 19981) for a sample of clouds, where it is considered that the pressure is produced by the surrounding 
diffuse ambient gas. Thus, if this clump is supported by external pressure it is required a pressure source, which in the present 
context it is likely to be the SNR shocks running into the outer layers of the molecular gas. The derived Pext Iks is v ery similar to 
the pr essure obtained towards a region associated with the face-on shock in the Vela SNR using FUSE observations (ISankrit et al.l 
l2Q0ll) . providing further support to our hypothesis that a SNR shock front contributes to maintain bound the observed ^^CO-clump. 

Finally we perform a SEP analysis t o the IR counterpart of Hii region G018.630-f0.30, the WISE sour ce J182348.08-123316 .1 
(ICutri & et al.ll2013l: IWright et~^l201Ql) . The fitting of the SED was done using the tool developed bv iRobitaille et al.l (l2007lPI 
considering the fluxes at the WISE 3.4, 4.6, 12, and 22 jjm bands, an Ai^ between 14 and 50 mag. (see Paper I) and a distance range 
of 13-15 kpc. The SED results show that the source is indeed a young massive star (age ~10^ yrs and mass ~ 12 M©), suggesting 
that its formation should be approximately coeval with the SN explosion. If the SED is performed using the near distance (a range 
1-3 kpc), we find that the source should have a lower mass (about 3 M©), which would be unable to ionize the gas and generate 
an Hii region. This is an indirect confirmation that the Hii region are located at the farther distance, suggesting that it should be 
embedded in the molecular cloud southwards the SNR. 


3.2. Region 2 

The emission from this region is significantly weaker than that detected towards Reg.l. By integrating the J=3-2 
emission between 15 and 30 km s“^ (see Eig.|4l-up), the whole velocity range where there is emission, we do not find any conspicuous 
morphological correspondence between the molecular gas and the most flattened border of the SNR which might suggests some 
kind of interaction. To appreciate in more detail the behavior of the molecular gas in this region, we present in Eig. 0-bottom the 
emission in a series of channel maps integrated in steps of 2 km s“^ It can be seen that the molecular gas has a very clumpy 
distribution not only in the plane of the sky but also along the line of sight. The emission shown in the three first panels in 
Eig.0bottom may correspond to gas located in front (along the line of sight) of the SNR G18.9-h0.3, while the emission displayed 
in the remaining panels would suggest a correspondence between the SNR and the molecular gas (see mainly the small clump 
towards the west at panel 24.5 and 26.5km s“^). We note that the molecular emission towards the bottom left corner of the surveyed 
region, mainly at 22.5 km s“^, should belong to the upper right border of the structure analysed in Reg. 1. Thus, even though the gas 

^ http ://caravan. astro.wisc.edu/protostars/ 
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distribution in the last panels of Fig.|4]-bottom may suggest a correspondence between the SNR and the molecular gas, these results 
allow us to conclude that the unusual flattened border of SNR G18.8+0.3 must have some physical origin other than interaction with 
dense environmental gas. 



Fig. 4. Left: J=3-2 emission integrated between 15 and 30 km s“^ with levels of 14, 18, and 28 K km s~^ (white contours) towards Reg.2. 

The color ranges are the same as in Fig.|2] Right: channel maps of the J=3-2 emission (red with white contours) integrated in steps of 2 km 
s“k The contours levels are 2.5, 4.0, and 6.0 K km s“^ 


4. Summary and concluding remarks 

We report the analysis of the molecular gas towards the southern environment of the SNR G18.8+0.3 using ASTE observations. We 
have surveyed two regions, one towards the SNR southernmost ‘corner’ which in addition harbors the Hit region GOl8.630+0.309, 
and a second one towards the most flattened border of the SNR. The main results can be summarized as follows: 

(a) In the first region, named Reg. 1/Reg. lb, we discovered a molecular feature with a good morphological correspondence with 
the SNR southernmost corner, where the indentations in the morphology of the SNR radio continuum emission are complemented by 
protrusions in the molecular cloud, strongly suggesting a SNR-molecular cloud interaction. Additionally, analysing this molecular 
cloud with the J=3-2 line we found a molecular feature composed by two clumps, not correlated with the peak, that 
partially surrounds the Hit region GOl8.630+0.309 and correlates with the thermal dust emission at 870 jum. We found that the 
virial mass of the ^^CO-clump considerably differ from the mass value obtained through other methods, suggesting that it is not 
gravitationally bound. Finally, from a SED fitting to the IR counterpart of the Hit region, we confirm that it is located at the far 
distance, suggesting that it should be embedded in the molecular cloud southwards the SNR. 

(b) The second region, named Reg.2, c overs a portion of th e most flattened border of SNR G18.8+0.3, where originally it was 
suggested an interaction with a dense cloud (iDubner et al.lll999h . Strikingly, our J=3-2 observations show a clumpy molecular 
structure without any morphological correspondence with the fiat border of the SNR that could suggest a causal relationship. The 
better angular resolution attained in these new data served to reveal that the fiat southern border of the SNR must originate in a 
different process other than the compression of a dense molecular cloud. 
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